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Floating Point Routines for the 6502 by Roy Rankin and Steve Wozniak 

Originally published in the August 1976 issue of Dr. Dobb's Journal, these 
floating point routines allow 6502 users to perform most of the more popular and 
desired floating point and transcendental functions, namely: Natural Log, Common 
Log, Addition, Subtraction, Multiplication, Division, and conversions between 
floating and fixed point numbers. 

Errata for Rankin's 6502 Floating Point Routines by Roy Rankin 

In the November/December issue of Dr. Dobb's Journal Roy Rankin published three 
error corrections to the Floating Point Routines presented above. 

Floating Point Implementation in the Apple II by Steve Wozniak 

An almost identical set of the above routines appeared in the original manual 
for the Apple II (the Red Book, January 1978) . Documentation for these routines 
appeared in another book, the Wozpak II, in November 1979. 



Floating Point Routines for the 6502 by Roy Rankin and Steve Wozniak 

Originally published in the August 1976 issue of Dr. Dobb's Journal, these 
floating point routines allow 6502 users to perform most of the more popular and 
desired floating point and transcendental functions, namely: Natural Log, Common 
Log, Addition, Subtraction, Multiplication, Division, and conversions between 
floating and fixed point numbers. 

Dr. Dobb's Journal, August 1976, pages 17-19. 

Floating Point Routines for the 6502 

by Roy Rankin, Department of Mechanical Engineering, 
Stanford University, Stanford, CA 94305 
(415) 497-1822 

and 

Steve Wozniak, Apple Computer Company 
770 Welch Road, Suite 154 
Palo Alto, CA 94304 
(415) 326-4248 

Editor's Note: Although these routines are for the 6502, it 
would appear that one could generate equivalent routines for 
most of the "traditional" microprocessors, relatively easily, 
by following the flow of the algorithms given in the excellent 
comments included in the program listing. This is particularly 
true of the transcendental functions, which were directly modeled 
after well-known and proven algorithms, and for which, the 
comments are relatively machine independent. 



These floating point routines allow 6502 users to perform 
most of the more popular and desired floating point and 
transcendental functions, namely: 

Natural Log - LOG 

Common Log - LOG10 

Exponential - EXP 

Floating Add - FADD 

Floating Subtract - FSUB 

Floating Multiply - FMUL 

Floating Divide - FDIV 

Convert Floating to Fixed - FIX 

Convert Fixed to Floating - FLOAT 

They presume a four-byte floating point operand consisting of 
a one-byte exponent ranging from -128 to +127 and a 
24-bit two's complement mantissa between 1.0 and 2.0. 

The floating point routines were done by Steve Wozniak, 
one of the principals in Apple Computer Company. The 
transcendental functions were patterned after those offered by 
Hewlett-Packard for their HP2100 minicomputer (with some 
modifications), and were done by Roy Rankin, a Ph.D. student 
at Stanford University. 

There are three error traps; two for overflow, and one for 
prohibited logarithm argument. ERROR (1D06) is the error 
exit used in the event of a non-positive log argument. OVFLW 
(1E3B) is the error exit for overflow occuring during calculation 
of e to some power. OVFL (1FE4) is the error exit for 
overflow in all of the floating point routines. There is no 
trap for underflow; in such cases, the result is set to 0.0. 

All routines are called and exited in a uniform manner: 

The arguments(s) are placed in the specified floating point 

storage locations (for specifics, see the documentation preceeding 

each routine in the listing), then a JSR is used to 

enter the desired routine. Upon normal completion, the 

called routine is exited via a subroutine return instruction (RTS) . 

Note: The preceeding documentation was written by the Editor, based 
on phone conversations with Roy and studying the listing. There is a 
high probability that it is correct. However, since it was not written 
nor reviewed by the authors of these routines, the preceeding 
documentation may contain errors in concept or in detail. 

-- JCW, Jr. 

In the Exponent: 
00 Represents -128 

7F Represents -1 

80 Represents 

81 Represents +1 

FF Represents +127 



Exponent Two's Complement Mantissa 
SEEEEEEE SM.IWIIWI MMIWIIWI MMIWIIWI 
n n+1 n+2 n+3 



* JULY 5, 1976 

* BASIC FLOATING POINT ROUTINES 

* FOR 6502 MICROPROCESSOR 

* BY R. RANKIN AND S. WOZNIAK 

* CONSISTING OF: 

* NATURAL LOG 

* COMMON LOG 

* EXPONENTIAL (E**X) 

* FLOAT FIX 

* FADD FSUB 

* FMUL FDIV 



FLOATING POINT REPRESENTATION (4-BYTES) 
EXPONENT BYTE 1 
MANTISSA BYTES 2-4 



0003 










ORG 


3 


0003 


EA 






SIGN 


NOP 




0004 


EA 






X2 


NOP 




0005 


00 


00 


00 


M2 


BSS 


3 


0008 


EA 






XI 


NOP 




0009 


00 


00 


00 


Ml 


BSS 


3 


000C 








E 


BSS 


4 


0010 








Z 


BSS 


4 


0014 








T 


BSS 


4 


0018 








SEXP 


BSS 


4 


001C 


00 






INT 

* 


BSS 


1 


1D00 








* 


ORG 


$1D00 










* 
* 
* 


NATURAL LOi 


1D00 


A5 


09 




LOG 


LDA 


Ml 


1D02 


F0 


02 






BEQ 


ERROR 



MANTISSA: TWO'S COMPLIMENT REPRESENTATION WITH SIGN IN 
MSB OF HIGH -ORDER BYTE. MANTISSA IS NORMALIZED WITH AN 
ASSUMED DECIMAL POINT BETWEEN BITS 5 AND 6 OF THE HIGH-ORDER 
BYTE. THUS THE MANTISSA IS IN THE RANGE 1. TO 2. EXCEPT 
WHEN THE NUMBER IS LESS THAN 2**(-128). 

EXPONENT: THE EXPONENT REPRESENTS POWERS OF TWO. THE 
REPRESENTATION IS 2'S COMPLIMENT EXCEPT THAT THE SIGN 
BIT (BIT 7) IS COMPLIMENTED. THIS ALLOWS DIRECT COMPARISON 
OF EXPONENTS FOR SIZE SINCE THEY ARE STORED IN INCREASING 
NUMERICAL SEQUENCE RANGING FROM $00 (-128) TO $FF (+127) 
($ MEANS NUMBER IS HEXADECIMAL). 

REPRESENTATION OF DECIMAL NUMBERS: THE PRESENT FLOATING 

POINT REPRESENTATION ALLOWS DECIMAL NUMBERS IN THE APPROXIMATE 
RANGE OF 10** (-38) THROUGH 10** (38) WITH 6 TO 7 SIGNIFICANT 
DIGITS. 



SET BASE PAGE ADRESSES 

EXPONENT 2 
MANTISSA 2 
EXPONENT 1 
MANTISSA 1 
SCRATCH 



STARTING LOCATION FOR LOG 
NATURAL LOG OF MANT/EXP1 WITH RESULT IN MANT/EXP1 



1D04 
1D06 

1D07 
1D0A 
1D0C 
1D0E 
1D10 
1D12 
1D14 
1D16 
1D18 
1D1B 
1D1D 
1D1F 
1D21 
1D23 
1D25 
1D28 
1D2A 
1D2B 
1D2D 
1D30 
1D32 
1D34 
1D36 
1D38 
1D3A 
1D3D 
1D3F 
1D40 
1D42 
1D45 
1D47 
1D49 
1D4B 
1D4C 
1D4E 
1D51 
1D53 
1D55 
1D57 
1D59 
1D5A 
1D5C 
1D5F 
1D62 
1D64 
1D67 
1D69 
1D6A 
1D6C 
1D6F 
1D71 
1D74 
1D76 
1D77 
1D79 
1D7C 
1D7E 
1D81 
1D83 
1D84 
1D86 



10 01 
00 

20 1C IF 
A5 04 
A0 80 

84 04 
49 80 

85 0A 
A9 00 
85 09 

20 2C IF 
A2 03 
B5 04 
95 10 
B5 08 
95 18 
BD Dl ID 
95 08 
CA 

10 F0 
20 4A IF 
A2 03 
B5 08 
95 14 
B5 10 
95 08 
BD Dl ID 
95 04 
CA 

10 F0 
20 50 IF 
A2 03 
B5 14 
95 04 
CA 

10 F9 
20 9D IF 
A2 03 
B5 08 
95 14 
95 04 
CA 

10 F7 
20 77 IF 
20 1C IF 
A2 03 
BD El ID 
95 08 
CA 

10 F8 
20 4A IF 
A2 03 
BD DD ID 
95 04 
CA 

10 F8 
20 9D IF 
A2 03 
BD D9 ID 
95 04 
CA 

10 F8 
20 50 IF 



ERROR 
CONT 



SEXP1 



SAVET 



TM2 



MIT 



MIC 



M2MB 



M2A1 



BPL CONT 
BRK 

JSR SWAP 

LDA X2 

LDY =$80 

STY X2 

EOR =$80 

STA Ml+1 

LDA =0 

STA Ml 

JSR FLOAT 

LDX =3 

LDA X2,X 

STA Z,X 

LDA XI, X 

STA SEXP,X 

LDA R22,X 

STA XI, X 

DEX 

BPL SEXP1 

JSR FSUB 

LDX =3 

LDA XI, X 

STA T,X 

LDA Z,X 

STA XI, X 

LDA R22,X 

STA X2,X 

DEX 

BPL SAVET 

JSR FADD 

LDX =3 

LDA T,X 

STA X2,X 

DEX 

BPL TM2 

JSR FDIV 

LDX =3 

LDA XI, X 

STA T,X 

STA X2,X 

DEX 

BPL MIT 

JSR FMUL 

JSR SWAP 

LDX =3 

LDA C,X 

STA XI, X 

DEX 

BPL MIC 

JSR FSUB 

LDX =3 

LDA MB,X 

STA X2,X 

DEX 

BPL M2MB 

JSR FDIV MB/(T*T-C) 

LDX =3 

LDA A1,X 

STA X2,X 

DEX 

BPL M2A1 

JSR FADD 



IF ARG>0 OK 
ERROR ARG<=0 

MOVE ARG TO EXP/MANT2 
HOLD EXPONENT 

SET EXPONENT 2 TO ($80) 

COMPLIMENT SIGN BIT OF ORIGINAL EXPONENT 

SET EXPONENT INTO MANTISSA 1 FOR FLOAT 

CLEAR MSB OF MANTISSA 1 
CONVERT TO FLOATING POINT 
4 BYTE TRANSFERS 

COPY MANTISSA TO Z 

SAVE EXPONENT IN SEXP 

LOAD EXP/MANT1 WITH SQRT(2) 



Z-SQRT(2) 

4 BYTE TRANSFER 

SAVE EXP/MANT1 AS T 

LOAD EXP/MANT1 WITH Z 

LOAD EXP/MANT2 WITH SQRT(2) 



Z+SQRT(2) 

4 BYTE TRANSFER 

LOAD T INTO EXP/MANT2 



T=(Z-SQRT(2))/(Z+SQRT(2)) 
4 BYTE TRANSFER 

COPY EXP/MANT1 TO T AND 
LOAD EXP/MANT2 WITH T 



T*T 

MOVE T*T TO EXP/MANT2 

4 BYTE TRANSFER 

LOAD EXP/MANT1 WITH C 



T*T-C 

4 BYTE TRANSFER 

LOAD EXP/MANT2 WITH MB 



LOAD EXP/MANT2 WITH Al 



MB/(T*T-C)+A1 



1D89 


A2 


03 






LDX 


=3 


4 BYTE TRANSFER 




1D8B 


B5 


14 




M2T 


LDA 


T,X 






1D8D 


95 


04 






STA 


X2,X 


LOAD EXP/MANT2 WITH T 




1D8F 


CA 








DEX 








1D90 


10 


F9 






BPL 


M2T 






1D92 


20 


77 


IF 




JSR 


FMUL 


(MB/(T*T-C)+A1)*T 




1D95 


A2 


03 






LDX 


=3 


4 BYTE TRANSFER 




1D97 


BD 


E5 


ID 


M2MHL 


LDA 


MHLF,X 






1D9A 


95 


04 






STA 


X2,X 


LOAD EXP/MANT2 WITH MHLF (.5) 




1D9C 


CA 








DEX 








1D9D 


10 


F8 






BPL 


M2MHL 






1D9F 


20 


50 


IF 




JSR 


FADD 


+ .5 




1DA2 


A2 


03 






LDX 


=3 


4 BYTE TRANSFER 




1DA4 


B5 


18 




LDEXP 


LDA 


SEXP,X 






1DA6 


95 


04 






STA 


X2,X 


LOAD EXP/MANT2 WITH ORIGINAL EXPOI* 


JEI 


1DA8 


CA 








DEX 








1DA9 


10 


F9 






BPL 


LDEXP 






1DAB 


20 


50 


IF 




JSR 


FADD 


+EXPN 




1DAE 


A2 


03 






LDX 


=3 


4 BYTE TRANSFER 




1DB0 


BD 


D5 


ID 


MLE2 


LDA 


LE2,X 






1DB3 


95 


04 






STA 


X2,X 


LOAD EXP/MANT2 WITH LN(2) 




1DB5 


CA 








DEX 








1DB6 


10 


F8 






BPL 


MLE2 






1DB8 


20 


77 


IF 




JSR 


FMUL 


*LN(2) 




1DBB 


60 






* 


RTS 




RETURN RESULT IN MANT/EXP1 












* 
* 


COMMON LOG OF MANT/EXP1 RESULT IN MANT/EXP1 




1DBC 


20 


00 


ID 


LOG10 


JSR 


LOG 


COMPUTE NATURAL LOG 




1DBF 


A2 


03 






LDX 


=3 






1DC1 


BD 


CD 


ID 


L10 


LDA 


LN10,X 






1DC4 


95 


04 






STA 


X2,X 


LOAD EXP/MANT2 WITH 1/LN(10) 




1DC6 


CA 








DEX 








1DC7 


10 


F8 






BPL 


L10 






1DC9 


20 


77 


IF 




JSR 


FMUL 


LOG10(X)=LN(X)/LN(10) 




1DCC 


60 






* 


RTS 








1DCD 


7E 


6F 




LN10 


DCM 


0.4342945 






2D 


ED 














1DD1 


80 


5A 




R22 


DCM 


1.4142136 SQRT(2) 






02 


7A 














1DD5 


7F 


58 




LE2 


DCM 


0.69314718 LOG BASE E OF 2 






B9 


0C 














1DD9 


80 


52 




Al 


DCM 


1.2920074 






80 


40 














1DDD 


81 


AB 




MB 


DCM 


-2.6398577 






86 


49 














1DE1 


80 


6A 




C 


DCM 


1.6567626 






08 


66 














1DE5 


7F 
00 


40 
00 




MHLF 

* 


DCM 


0.5 






1E00 








* 


ORG 


$1E00 


STARTING LOCATION FOR EXP 












* 
* 


EXP OF MANT/EXP1 RESULT IN MANT/EXP1 




1E00 


A2 


03 




EXP 


LDX 


=3 


4 BYTE TRANSFER 




1E02 


BD 


D8 


IE 




LDA 


L2E,X 






1E05 


95 


04 






STA 


X2,X 


LOAD EXP/MANT2 WITH LOG BASE 2 OF 


E 


1E07 


CA 








DEX 








1E08 


10 


F8 






BPL 


EXP+2 






1E0A 


20 


77 


IF 




JSR 


FMUL 


L0G2(3)*X 




1E0D 


A2 


03 






LDX 


=3 


4 BYTE TRANSFER 




1E0F 


B5 


08 




FSA 


LDA 


XI, X 







1E11 
1E13 
1E14 
1E16 
1E19 
1E1B 
1E1D 
1E1E 
1E20 
1E22 
1E24 
1E26 
1E27 
1E29 
1E2B 
1E2D 
1E2F 
1E31 
1E33 
1E35 
1E37 
1E38 
1E3A 



1E3C 
1E3F 
1E41 
1E43 
1E45 
1E46 
1E48 
1E4B 
1E4D 
1E4F 
1E51 
1E53 
1E54 
1E56 
1E59 
1E5B 
1E5E 
1E60 
1E62 
1E64 
1E65 
1E67 
1E6A 
1E6C 
1E6F 
1E71 
1E72 
1E74 
1E77 
1E79 
1E7B 
1E7D 
1E80 
1E82 
1E84 
1E86 
1E87 
1E89 



95 10 
CA 

10 F9 
20 E8 IF 
A5 0A 
85 1C 
38 

E9 7C 
A5 09 
E9 00 
10 15 
18 

A5 0A 
69 78 
A5 09 
69 00 
10 0B 
A9 00 
A2 03 
95 08 
CA 

10 FB 
60 



1E3B 00 



20 2C IF 
A2 03 
B5 10 
95 04 
CA 

10 F9 
20 4A IF 
A2 03 
B5 08 
95 10 
95 04 
CA 

10 F7 
20 77 IF 
A2 03 
BD DC IE 
95 04 
B5 08 
95 18 
CA 

10 F4 
20 50 IF 
A2 03 
BD E0 IE 
95 04 
CA 

10 F8 
20 9D IF 
A2 03 
B5 08 
95 14 
BD E4 IE 
95 08 
B5 18 
95 04 
CA 

10 F0 
20 77 IF 



ZERO 



* 
OVFLW 



STA Z,X 

DEX 

BPL FSA 

JSR FIX 

LDA Ml+1 

STA INT 

SEC 

SBC =124 

LDA Ml 

SBC =0 

BPL OVFLW 

CLC 

LDA Ml+1 

ADC =120 

LDA Ml 

ADC =0 

BPL CONTIN 

LDA =0 

LDX =3 

STA XI, X 

DEX 

BPL ZERO 

RTS 

BRK 



CONTIN JSR FLOAT 
ENTD 



LDX 
LDA 



=3 
Z,X 



ZSAV 



LA2 



LB2 



DLOAD 



STA X2,X 
DEX 

BPL ENTD 
JSR FSUB 
LDX =3 
LDA XI, X 
STA Z,X 
STA X2,X 
DEX 

BPL ZSAV 
JSR FMUL 
LDX =3 
LDA A2,X 
STA X2,X 
LDA XI, X 
STA SEXP,X 
DEX 

BPL LA2 
JSR FADD 
LDX =3 
LDA B2,X 
STA X2,X 
DEX 

BPL LB2 
JSR FDIV 
LDX =3 
LDA XI, X 
STA T,X 
LDA C2,X 
STA XI, X 
LDA SEXP,X 
STA X2,X 
DEX 

BPL DLOAD 
JSR FMUL 



STORE EXP/MANT1 IN Z 

SAVE Z=LN(2)*X 

CONVERT CONTENTS OF EXP/MANT1 TO AN INTEGER 

SAVE RESULT AS INT 

SET CARRY FOR SUBTRACTION 

INT- 124 



OVERFLOW INT>=124 
CLEAR CARRY FOR ADD 

ADD 120 TO INT 



IF RESULT POSITIVE CONTINUE 

INT<-120 SET RESULT TO ZERO AND RETURN 

4 BYTE MOVE 

SET EXP/MANT1 TO ZERO 



RETURN 
OVERFLOW 
FLOAT INT 

LOAD EXP/MANT2 WITH Z 



Z*Z-FLOAT(INT) 
4 BYTE MOVE 

SAVE EXP/MANT1 IN Z 

COPY EXP/MANT1 TO EXP/MANT2 



Z*Z 

4 BYTE MOVE 

LOAD EXP/MANT2 WITH A2 

SAVE EXP/MANT1 AS SEXP 



Z*Z+A2 

4 BYTE MOVE 

LOAD EXP/MANT2 WITH B2 



T=B/(Z*Z+A2) 
4 BYTE MOVE 

SAVE EXP/MANT1 AS T 

LOAD EXP/MANT1 WITH C2 

LOAD EXP/MANT2 WITH SEXP 

Z*Z*C2 



1E8C 
1E8F 
1E91 
1E93 
1E95 
1E96 
1E98 
1E9B 
1E9D 
1EA0 
1EA2 
1EA3 
1EA5 
1EA8 
1EAB 
1EAD 
1EAF 
1EB1 
1EB2 
1EB4 
1EB7 
1EB9 
1EBB 
1EBD 
1EBE 
1EC0 
1EC3 
1EC5 
1EC8 
1ECA 
1ECB 
1ECD 
1ED0 
1ED1 
1ED3 
1ED5 
1ED7 
1ED8 

1EDC 

1EE0 

1EE4 

1EE8 



20 1C IF 
A2 03 
B5 14 
95 08 
CA 

10 F9 
20 4A IF 
A2 03 
BD E8 IE 
95 04 
CA 

10 F8 
20 50 IF 
20 1C IF 
A2 03 
B5 10 
95 08 
CA 

10 F9 
20 4A IF 
A2 03 
B5 10 
95 04 
CA 

10 F9 
20 9D IF 
A2 03 
BD E5 ID 
95 04 
CA 

10 F8 
20 50 IF 
38 

A5 1C 
65 08 

85 08 
60 

80 5C 
55 IE 

86 57 
6A El 
89 4D 
3F ID 
7B 46 
FA 70 
83 4F 
A3 03 



JSR 
LDX 

LTMP LDA 
STA 
DEX 
BPL 
JSR 
LDX 

LDD LDA 
STA 
DEX 
BPL 
JSR 
JSR 
LDX 

LFA LDA 
STA 
DEX 
BPL 
JSR 
LDX 

LF3 LDA 
STA 
DEX 
BPL 
JSR 
LDX 

LD12 LDA 
STA 
DEX 
BPL 
JSR 
SEC 
LDA 
ADC 
STA 
RTS 

L2E DCM 



SWAP 
=3 
T,X 
XI, X 

LTMP 

FSUB 

=3 

D,X 

X2,X 

LDD 

FADD 

SWAP 

=3 

Z,X 

XI, X 

LFA 

FSUB 

=3 

Z,X 

X2,X 

LF3 

FDIV 

=3 

MHLF,X 

X2,X 

LD12 
FADD 



MOVE EXP/MANT1 TO EXP/MANT2 
4 BYTE TRANSFER 

LOAD EXP/MANT1 WITH T 



C2*Z*Z-B2/(Z*Z+A2) 
4 BYTE TRANSFER 

LOAD EXP/MANT2 WITH D 



D+C2*Z*Z-B2/(Z*Z+A2) 

MOVE EXP/MANT1 TO EXP/MANT2 

4 BYTE TRANSFER 

LOAD EXP/MANT1 WITH Z 



- Z+D+C2*Z*Z- B2/ ( Z*Z+A2 ) 
4 BYTE TRANSFER 

LOAD EXP/MANT2 WITH Z 



2/(**** ) 

4 BYTE TRANSFER 

LOAD EXP/MANT2 WITH .5 



A2 
B2 
C2 
D 



DCM 
DCM 
DCM 
DCM 



+Z/(***)+.5 

ADD INT TO EXPONENT WITH CARRY SET 
INT TO MULTIPLY BY 
XI 2**(INT+1) 
XI RETURN RESULT TO EXPONENT 

RETURN ANS= ( . 5+Z/ ( - Z+D+C2*Z*Z- B2/ ( Z*Z+A2 ) ) *2** ( INT+1 ) 
1.4426950409 LOG BASE 2 OF E 

87.417497202 

617.9722695 

.03465735903 

9.9545957821 



BASIC FLOATING POINT ROUTINES 



1F00 
1F00 
1F01 
1F03 
1F05 
1F07 
1F09 
1F0A 
1F0C 
1F0D 
1F0F 
1F12 
1F14 



18 

A2 02 
B5 09 
75 05 
95 09 
CA 

10 F7 
60 

06 03 
20 12 IF 
24 09 
10 05 



ADD 
ADD1 



ORG $1F00 

CLC 

LDX =$02 

LDA M1,X 

ADC M2,X 

STA M1,X 

DEX 

BPL ADD1 

RTS 

ASL SIGN 



START OF BASIC FLOATING POINT ROUTINES 

CLEAR CARRY 

INDEX FOR 3- BYTE ADD 

ADD A BYTE OF MANT2 TO MANT1 



ADVANCE INDEX TO NEXT MORE SIGNIF.BYTE 
LOOP UNTIL DONE. 
RETURN 
MD1 ASL SIGN CLEAR LSB OF SIGN 

JSR ABSWAP ABS VAL OF MANT1, THEN SWAP MANT2 
ABSWAP BIT Ml MANT1 NEG? 

BPL ABSWP1 NO, SWAP WITH MANT2 AND RETURN 



1F16 20 8F IF JSR FCOMPL YES, COMPLIMENT IT. 

1F19 E6 03 INC SIGN INCR SIGN, COMPLEMENTING LSB 

1F1B 38 ABSWP1 SEC SET CARRY FOR RETURN TO MUL/DIV 



SWAP EXP/MANT1 WITH EXP/MANT2 



1F1C A2 04 SWAP LDX =$04 INDEX FOR 4- BYTE SWAP. 

1F1E 94 0B SWAP1 STY E-l.X 

1F20 B5 07 LDA Xl-l.X SWAP A BYTE OF EXP/MANT1 WITH 

1F22 B4 03 LDY X2-1.X EXP/MANT2 AND LEAVEA COPY OF 

1F24 94 07 STY Xl-l.X MANT1 IN E(3BYTES). E+3 USED. 

1F26 95 03 STA X2-1.X 

1F28 CA DEX ADVANCE INDEX TO NEXT BYTE 

1F29 DO F3 BNE SWAP1 LOOP UNTIL DONE. 

1F2B 60 RTS 

* CONVERT 16 BIT INTEGER IN Ml(HIGH) AND Ml+KLOW) TO F.P. 

* RESULT IN EXP/MANT1. EXP/MANT2 UNEFFECTED 

1F2C A9 8E FLOAT LDA =$8E 

1F2E 85 08 STA XI SET EXPN TO 14 DEC 

1F30 A9 00 LDA =0 CLEAR LOW ORDER BYTE 

1F32 85 0B STA Ml+2 

1F34 F0 08 BEQ NORM NORMALIZE RESULT 

1F36 C6 08 N0RM1 DEC XI DECREMENT EXP1 

1F38 06 0B ASL Ml+2 

1F3A 26 0A ROL Ml+1 SHIFT MANT1 (3 BYTES) LEFT 

1F3C 26 09 ROL Ml 

1F3E A5 09 NORM LDA Ml HIGH ORDER MANT1 BYTE 

1F40 0A ASL UPPER TWO BITS UNEQUAL? 

1F41 45 09 EOR Ml 

1F43 30 04 BMI RTS1 YES, RETURN WITH MANT1 NORMALIZED 

1F45 A5 08 LDA XI EXP1 ZERO? 

1F47 DO ED BNE N0RM1 NO, CONTINUE NORMALIZING 

1F49 60 RTS1 RTS RETURN 

* EXP/MANT2-EXP/MANT1 RESULT IN EXP/MANT1 

1F4A 20 8F IF FSUB JSR FCOMPL CMPL MANT1 CLEARS CARRY UNLESS ZERO 

1F4D 20 5D IF SWPALG JSR ALGNSW RIGHT SHIFT MANT1 OR SWAP WITH MANT2 ON CARRY 

* ADD EXP/MANT1 AND EXP/MANT2 RESULT IN EXP/MANT1 

1F50 A5 04 FADD LDA X2 

1F52 C5 08 CMP XI COMPARE EXP1 WITH EXP2 

1F54 DO F7 BNE SWPALG IF UNEQUAL, SWAP ADDENDS OR ALIGN MANTISSAS 

1F56 20 00 IF JSR ADD ADD ALIGNED MANTISSAS 

1F59 50 E3 ADDEND BVC NORM NO OVERFLOW, NORMALIZE RESULTS 

1F5B 70 05 BVS RTLOG OV: SHIFT MANT1 RIGHT. NOTE CARRY IS CORRECT SIGN 

1F5D 90 BD ALGNSW BCC SWAP SWAP IF CARRY CLEAR, ELSE SHIFT RIGHT ARITH. 

1F5F A5 09 RTAR LDA Ml SIGN OF MANT1 INTO CARRY FOR 

1F61 0A ASL RIGHT ARITH SHIFT 

1F62 E6 08 RTLOG INC XI INCR EXP1 TO COMPENSATE FOR RT SHIFT 

1F64 F0 7E BEQ OVFL EXP1 OUT OF RANGE. 

1F66 A2 FA RTL0G1 LDX =$FA INDEX FOR 6 BYTE RIGHT SHIFT 

1F68 A9 80 R0R1 LDA =$80 

1F6A B0 01 BCS R0R2 

1F6C 0A ASL 

1F6D 56 0F R0R2 LSR E+3,X SIMULATE ROR E+3,X 

1F6F 15 0F ORA E+3,X 



1F71 95 OF 

1F73 E8 

1F74 DO F2 

1F76 60 



1F77 
1F7A 
1F7C 
1F7F 
1F80 
1F83 
1F85 
1F88 
1F89 
1F8B 
1F8D 
1F8F 
1F90 
1F92 
1F94 
1F96 
1F98 
1F99 
1F9B 



20 0D 
65 08 
20 CD 
18 

20 66 
90 03 
20 00 
88 

10 F5 
46 03 
90 AF 
38 

A2 03 
A9 00 
F5 08 
95 08 
CA 

D0 F7 
F0 BC 



IF 



IF 



IF 



IF 



STA E+3,X 

INX 

BNE ROR1 

RTS 



NEXT BYTE OF SHIFT 
LOOP UNTIL DONE 
RETURN 



* 

FMUL JSR 
ADC 
JSR 
CLC 

MUL1 JSR 
BCC 
JSR 

MUL2 DEY 
BPL 

MDEND LSR 

NORMX BCC 

FCOMPL SEC 
LDX 

COMPL1 LDA 
SBC 
STA 
DEX 
BNE 
BEQ 



EXP/MANT1 X EXP/MANT2 RESULT IN EXP/MANT1 



MD1 ABS. VAL OF MANT1, MANT2 

XI ADD EXP1 TO EXP2 FOR PRODUCT EXPONENT 

MD2 CHECK PRODUCT EXP AND PREPARE FOR MUL 

CLEAR CARRY 
RTLOG1 MANT1 AND E RIGHT. (PRODUCT AND MPLIER) 
MUL2 IF CARRY CLEAR, SKIP PARTIAL PRODUCT 
ADD ADD MULTIPLICAN TO PRODUCT 

NEXT MUL ITERATION 
MUL1 LOOP UNTIL DONE 
SIGN TEST SIGN (EVEN/ODD) 
NORM IF EXEN, NORMALIZE PRODUCT, ELSE COMPLEMENT 

SET CARRY FOR SUBTRACT 
=$03 INDEX FOR 3 BYTE SUBTRACTION 
=$00 CLEAR A 
XI, X SUBTRACT BYTE OF EXP1 
XI, X RESTORE IT 

NEXT MORE SIGNIFICANT BYTE 
COMPL1 LOOP UNTIL DONE 
ADDEND NORMALIZE (OR SHIFT RIGHT IF OVERFLOW) 



1F9D 
1FA0 
1FA2 
1FA5 
1FA6 
1FA8 
1FAA 
1FAC 
1FAD 
1FAE 
1FB0 
1FB2 
1FB3 
1FB5 
1FB7 
1FB8 
1FBA 
1FBC 
1FBE 
1FC0 
1FC2 
1FC4 
1FC6 
1FC8 
1FC9 
1FCB 
1FCD 
1FCF 
1FD1 
1FD3 
1FD5 
1FD7 
1FD8 



20 OD 
E5 08 
20 CD 
38 

A2 02 
B5 05 
F5 OC 
48 
CA 

10 F8 
A2 FD 
68 

90 02 
95 08 
E8 

DO F8 
26 OB 
26 OA 
26 09 
06 07 
26 06 
26 05 
BO 1C 
88 

DO DA 
FO BE 
86 OB 
86 OA 
86 09 
BO OD 
30 04 
68 
68 



* 
IF FDIV 



EXP/MANT2 / EXP/MANT1 RESULT IN EXP/MANT1 



IF 



DIV1 



DIV2 



DIV3 



DIV4 



MD2 



JSR 
SBC 
JSR 
SEC 
LDX 
LDA 
SBC 
PHA 
DEX 
BPL 
LDX 
PLA 
BCC 
STA 
INX 
BNE 
ROL 
ROL 
ROL 
ASL 
ROL 
ROL 
BCS 
DEY 
BNE 
BEQ 
STX 
STX 
STX 
BCS 
BMI 
PLA 
PLA 



MD1 

XI 

MD2 

=$02 
M2,X 
E,X 



DIV2 
=$FD 

DIV4 
M2+3,X 

DIV3 

Ml+2 

Ml+1 

Ml 

M2+2 

M2+1 

M2 

OVFL 

DIV1 

MDEND 

Ml+2 

Ml+1 

Ml 

OVCHK 

MD3 



TAKE ABS VAL OF MANT1, MANT2 
SUBTRACT EXP1 FROM EXP2 
SAVE AS QUOTIENT EXP 
SET CARRY FOR SUBTRACT 
INDEX FOR 3- BYTE INSTRUCTION 

SUBTRACT A BYTE OF E FROM MANT2 

SAVE ON STACK 

NEXT MORE SIGNIF BYTE 

LOOP UNTIL DONE 

INDEX FOR 3- BYTE CONDITIONAL MOVE 

PULL A BYTE OF DIFFERENCE OFF STACK 

IF MANT2<E THEN DONT RESTORE MANT2 

NEXT LESS SIGNIF BYTE 
LOOP UNTIL DONE 

ROLL QUOTIENT LEFT, CARRY INTO LSB 



SHIFT DIVIDEND LEFT 

OVERFLOW IS DUE TO UNNORMALIZED DIVISOR 
NEXT DIVIDE ITERATION 
LOOP UNTIL DONE 23 ITERATIONS 
NORMALIZE QUOTIENT AND CORRECT SIGN 

CLR MANT1 (3 BYTES) FOR MUL/DIV 

IF EXP CALC SET CARRY, CHECK FOR OVFL 
IF NEG NO UNDERFLOW 
POP ONE 
RETURN LEVEL 



1FD9 


90 


B2 




BCC NORMX 


1FDB 


49 


80 


MD3 


EOR =$80 


1FDD 


85 


08 




STA XI 


1FDF 


A0 


17 




LDY =$17 


1FE1 


60 






RTS 


1FE2 


10 


F7 


OVCHK 


BPL MD3 


1FE4 


00 




OVFL 

* 
* 
* 
* 


BRK 

CONVERT EXI 
EXP/MANT2 


1FE5 


20 


5F IF 




JSR RTAR 


1FE8 


A5 


08 


FIX 


LDA XI 


1FEA 


C9 


8E 




CMP =$8E 


1FEC 


D0 


F7 




BNE FIX-3 


1FEE 


60 




RTRN 


RTS 

END 



CLEAR XI AND RETURN 

COMPLIMENT SIGN BIT OF EXP 

STORE IT 

COUNT FOR 24 MUL OR 23 DIV ITERATIONS 

RETURN 

IF POS EXP THEN NO OVERFLOW 



(HIGH) AND Ml+KLOW) 



SHIFT MANT1 RT AND INCREMENT EXPNT 

CHECK EXPONENT 

IS EXPONENT 14? 

NO, SHIFT 

RETURN 



Errata for Rankin's 6502 Floating Point Routines by Roy Rankin 

In the November/December issue of Dr. Dobb's Journal Roy Rankin published three 
error corrections to the Floating Point Routines presented above. 

Dr. Dobb's Journal, November/December 1976, page 57. 

ERRATA FOR RANKIN'S 6502 
FLOATING POINT ROUTINES 

Sept. 22, 1976 

Dear Jim, 

Subsequent to the publication of "Floating Point 
Routines for the 6502" (Vol.1, No. 7) an error which I made in 
the LOG routine came to light which causes improper results 
if the argument is less than 1. The following changes will 
correct the error. 



After: 






CONT 


JSR SWAP 


Add: 


A2 


00 




LDX =0 


After: 








STA Ml+1 


Delete: 








LDA =0 

STA Ml 


Add: 


10 
CA 


01 




BPL *+3 
DEX 




86 


09 




STX Ml 



;ido7) 
load x for high byte of exponent 

;iD12) 



IS EXPONENT NEGATIVE 

YES, SET X TO $FF 

SET UPPER BYTE OF EXPONENT 



3. Changes 1 and 2 shift the code by 3 bytes so add 3 to the 

addresses of the constants LN10 through MHLF whenever 

they are referenced. For example the address of LN10 changes 

from 1DCD to 1DD0. Note also that the entry point for 

LOG10 becomes 1DBF. The routines stays within the page 

and hence the following routines (EXP etc.) are not affected. 



Yours truly, 



Roy Rankin 

Dep. of Mech. Eng. 

Stanford University 



Floating Point Implementation in the Apple II by Steve Wozniak 

An almost identical set of the above routines appeared in the original manual 
for the Apple II (the Red Book, January 1978) . Documentation for these routines 
appeared in another book, the Wozpak II, in November 1979. 

Woz 6502 Floating Point Routines 

Apple II Reference Manual (Red Book), January 1978, pages 94-95. 

* APPLE- II FLOATING * 

* POINT ROUTINES * 

* * 

* COPYRIGHT 1977 BY * 

* APPLE COMPUTER INC. * 



* ALL RIGHTS RESERVED * 



F425 
F426 
F428 
F42A 
F42C 
F42E 
F42F 
F431 
F432 
F434 
F437 
F439 
F43B 
F43E 
F440 
F441 
F443 
F445 
F447 
F449 
F44B 
F44D 



F4 



18 

A2 02 
B5 F9 
75 F5 
95 F9 
CA 

10 F7 
60 

06 F3 
20 37 
24 F9 
10 05 
20 A4 F4 
E6 F3 
38 

A2 04 
94 FB 
B5 F7 
B4 F3 

94 F7 

95 F3 
CA 



S. WOZNIAK 






TITLE 
SIGN 
X2 

M2 

XI 

Ml 

E 

OVLOC 

ADD 

ADD1 



"FLOATING POINT 
EPZ $F3 



ROUTINES" 



MD1 
ABSWAP 



ABSWAP1 

SWAP 

SWAP1 



EPZ 
EPZ 
EPZ 
EPZ 
EPZ 
EQU 
ORG 
CLC 
LDX 
LDA 
ADC 
STA 
DEX 
BPL 
RTS 
ASL 
JSR 
BIT 
BPL 
JSR 
INC 
SEC 
LDX 
STY 
LDA 
LDY 
STY 
STA 
DEX 



$F4 

$F5 

$F8 

$F9 

$FC 

$3F5 

$F425 

#$2 
M1,X 
M2,X 
M1,X 

ADD1 

SIGN 

ABSWAP 

Ml 

ABSWAP1 

FCOMPL 

SIGN 

#$4 

E-1,X 

Xl-l.X 

X2-1.X 

Xl-l.X 

X2-1,X 



CLEAR CARRY 

INDEX FOR 3 -BYTE ADD. 

ADD A BYTE OF MANT2 TO MANT1 

INDEX TO NEXT MORE SIGNIF. BYTE. 

LOOP UNTIL DONE. 

RETURN 

CLEAR LSB OF SIGN. 

ABS VAL OF Ml, THEN SWAP WITH M2 

MANT1 NEGATIVE? 

NO, SWAP WITH MANT2 AND RETURN. 

YES, COMPLEMENT IT. 

INCR SIGN, COMPLEMENTING LSB. 

SET CARRY FOR RETURN TO MUL/DIV. 

INDEX FOR 4 BYTE SWAP. 

SWAP A BYTE OF EXP/MANT1 WITH 
EXP/MANT2 AND LEAVE A COPY OF 
MANT1 IN E (3 BYTES). E+3 USED 

ADVANCE INDEX TO NEXT BYTE 



F44E: 
F450: 
F451: 
F453: 
F455: 
F457: 
F459: 
F45B: 
F45D: 
F45F: 
F461: 
F463: 
F465: 
F467: 
F468: 
F46B: 
F46E: 
F470: 
F472: 
F474: 
F477: 
F479: 
F47B: 

F47D: 
F47F: 
F480: 
F482: 
F484: 
F486: 
F488: 
F489: 
F48B: 
F48C: 
F48F: 
F491: 
F494: 
F495: 
F498: 
F49A: 
F49D: 
F49E: 
F4A0: 
F4A2: 
F4A4: 
F4A5: 
F4A7: 
F4A9: 
F4AB: 
F4AD: 
F4AE: 
F4B0: 
F4B2: 
F4B5: 
F4B7: 
F4BA: 
F4BB: 
F4BD: 
F4BF: 
F4C1: 
F4C2: 
F4C3: 
F4C5: 
F4C7: 



F4 



D0 F3 
60 

A9 8E 
85 F8 
A5 F9 
C9 C0 
30 0C 
C6 F8 
06 FB 
26 FA 
26 F9 
A5 F8 
D0 EE 
60 

20 A4 F4 
20 7B F4 
A5 F4 
C5 F8 
D0 F7 
20 25 
50 EA 
70 05 
90 C4 

A5 F9 

0A 

E6 F8 

F0 75 

A2 FA 

76 FF 

E8 

D0 FB 

60 

20 32 

65 F8 

20 E2 

18 

20 84 F4 

90 03 

20 25 



FLOAT 

NORM1 



F4 



F4 



88 

10 F5 
46 F3 
90 BF 
38 

A2 03 
A9 00 
F5 F8 
95 F8 
CA 

D0 F7 
F0 C5 
20 32 
E5 F8 
20 E2 
38 

A2 02 
B5 F5 
F5 FC 
48 
CA 

10 F8 
A2 FD 
68 



F4 



NORM 

RTS1 
FSUB 
SWPALGN 
FADD 



ADDEND 

ALGNSWP 

RTAR 

RTLOG 

RTLOG1 
ROR1 



F4 FMUL 



MUL1 



MUL2 

MDEND 
NORMX 
FCOMPL 

COMPL1 



F4 FDIV 



DIV1 
DIV2 



DIV3 



BNE 
RTS 
LDA 
STA 
LDA 
CMP 
BMI 
DEC 
ASL 
ROL 
ROL 
LDA 
BNE 
RTS 
JSR 
JSR 
LDA 
CMP 
BNE 
JSR 
BVC 
BVS 
BCC 
ELSE SH 
LDA 
ASL 
INC 
BEQ 
LDX 
ROR 
I NX 
BNE 
RTS 
JSR 
ADC 
JSR 
CLC 
JSR 
BCC 
JSR 
DEY 
BPL 
LSR 
BCC 
SEC 
LDX 
LDA 
SBC 
STA 
DEX 
BNE 
BEQ 
JSR 
SBC 
JSR 
SEC 
LDX 
LDA 
SBC 
PHA 
DEX 
BPL 
LDX 
PLA 



SWAP1 

#$8E 

XI 

Ml 

#$C0 

RTS1 

XI 

Ml+2 

Ml+1 

Ml 

XI 

NORM1 

FCOMPL 

ALGNSWP 

X2 

XI 

SWPALGN 

ADD 

NORM 

RTLOG 

SWAP 

IFT RIGHT ARITH. 



LOOP UNTIL DONE. 

RETURN 

INIT EXP1 TO 14, 

THEN NORMALIZE TO FLOAT. 

HIGH-ORDER MANT1 BYTE. 

UPPER TWO BITS UNEQUAL? 

YES, RETURN WITH MANT1 NORMALIZED 

DECREMENT EXP1. 

SHIFT MANT1 (3 BYTES) LEFT. 

EXP1 ZERO? 

NO, CONTINUE NORMALIZING. 

RETURN. 

CMPL MANT1, CLEARS CARRY UNLESS 

RIGHT SHIFT MANT1 OR SWAP WITH 

COMPARE EXP1 WITH EXP2. 

IF #,SWAP ADDENDS OR ALIGN MANTS. 

ADD ALIGNED MANTISSAS. 

NO OVERFLOW, NORMALIZE RESULT. 

OV: SHIFT Ml RIGHT, CARRY INTO SIGN 

SWAP IF CARRY CLEAR, 



Ml SIGN OF MANT1 INTO CARRY FOR 

RIGHT ARITH SHIFT. 
XI INCR XI TO ADJUST FOR RIGHT SHIFT 
OVFL EXP1 OUT OF RANGE. 
#$FA INDEX FOR 6: BYTE RIGHT SHIFT. 
E+3,X 

NEXT BYTE OF SHIFT. 
ROR1 LOOP UNTIL DONE. 

RETURN. 
MD1 ABS VAL OF MANT1, MANT2 
XI ADD EXP1 TO EXP2 FOR PRODUCT EXP 
MD2 CHECK PROD. EXP AND PREP. FOR MUL 

CLEAR CARRY FOR FIRST BIT. 
RTLOG1 Ml AND E RIGHT (PROD AND MPLIER) 
MUL2 IF CARRY CLEAR, SKIP PARTIAL PROD 
ADD ADD MULTIPLICAND TO PRODUCT. 

NEXT MUL ITERATION. 
MUL1 LOOP UNTIL DONE. 
SIGN TEST SIGN LSB. 
NORM IF EVEN, NORMALIZE PROD, ELSE COMP 

SET CARRY FOR SUBTRACT. 
#$3 INDEX FOR 3 BYTE SUBTRACT. 
#$0 CLEAR A. 
XI, X SUBTRACT BYTE OF EXP1 . 
XI, X RESTORE IT. 

NEXT MORE SIGNIFICANT BYTE. 
COMPL1 LOOP UNTIL DONE. 
ADDEND NORMALIZE (OR SHIFT RT IF OVFL). 
MD1 TAKE ABS VAL OF MANT1, MANT2 . 
XI SUBTRACT EXP1 FROM EXP2 . 
MD2 SAVE AS QUOTIENT EXP. 

SET CARRY FOR SUBTRACT. 
#$2 INDEX FOR 3 -BYTE SUBTRACTION. 
M2,X 
E,X SUBTRACT A BYTE OF E FROM MANT2 . 

SAVE ON STACK. 

NEXT MORE SIGNIFICANT BYTE. 
DIV2 LOOP UNTIL DONE. 
#$FD INDEX FOR 3-BYTE CONDITIONAL MOVE 

PULL BYTE OF DIFFERENCE OFF STACK 



F4C8 


90 


02 




BCC 


DIV4 


F4CA 


95 


F8 




STA 


M2+3,X 


F4CC 


E8 




DIV4 


I NX 




F4CD 


D0 


F8 




BNE 


DIV3 


F4CF 


26 


FB 




ROL 


Ml+2 


F4D1 


26 


FA 




ROL 


Ml+1 


F4D3 


26 


F9 




ROL 


Ml 


F4D5 


06 


F7 




ASL 


M2+2 


F4D7 


26 


F6 




ROL 


M2+1 


F4D9 


26 


F5 




ROL 


M2 


F4DB 


B0 


1C 




BCS 


OVFL 


F4DD 


88 






DEY 




F4DE 


D0 


DA 




BNE 


DIV1 


F4E0 


F0 


BE 




BEQ 


MDEND 


F4E2 


86 


FB 


MD2 


STX 


Ml+2 


F4E4 


86 


FA 




STX 


Ml+1 


F4E6 


86 


F9 




STX 


Ml 


F4E8 


B0 


0D 




BCS 


OVCHK 


F4EA 


30 


04 




BMI 


MD3 


F4EC 


68 






PLA 




F4ED 


68 






PLA 




F4EE 


90 


B2 




BCC 


NORMX 


F4F0 


49 


80 


MD3 


EOR 


#$80 


F4F2 


85 


F8 




STA 


XI 


F4F4 


A0 


17 




LDY 


#$17 


F4F6 


60 






RTS 




F4F7 


10 


F7 


OVCHK 


BPL 


MD3 


F4F9 


4C 


F5 03 


OVFL 


JMP 
ORG 


OVLOC 
$F63D 


F63D 


20 


7D F4 


FIX1 


JSR 


RTAR 


F640 


A5 


F8 


FIX 


LDA 


XI 


F642 


10 


13 




BPL 


UNDFL 


F644 


C9 


8E 




CMP 


#$8E 


F646 


D0 


F5 




BNE 


FIX1 


F648 


24 


F9 




BIT 


Ml 


F64A 


10 


0A 




BPL 


FIXRTS 


F64C 


A5 


FB 




LDA 


Ml+2 


F64E 


F0 


06 




BEQ 


FIXRTS 


F650 


E6 


FA 




INC 


Ml+1 


F652 


D0 


02 




BNE 


FIXRTS 


F654 


E6 


F9 




INC 


Ml 


F656 


60 




FIXRTS 


RTS 




F657 


A9 


00 


UNDFL 


LDA 


#$0 


F659 


85 


F9 




STA 


Ml 


F65B 


85 


FA 




STA 


Ml+1 


F65D 


60 






RTS 





IF M2<E THEN DON'T RESTORE M2 . 

NEXT LESS SIGNIFICANT BYTE. 
LOOP UNTIL DONE. 

ROLL QUOTIENT LEFT, CARRY INTO LSB 



SHIFT DIVIDEND LEFT 

OVFL IS DUE TO UNNORMED DIVISOR 
NEXT DIVIDE ITERATION. 
LOOP UNTIL DONE 23 ITERATIONS. 
NORM. QUOTIENT AND CORRECT SIGN. 

CLEAR MANT1 (3 BYTES) FOR MUL/DIV. 

IF CALC. SET CARRY, CHECK FOR OVFL 
IF NEG THEN NO UNDERFLOW. 
POP ONE RETURN LEVEL. 

CLEAR XI AND RETURN. 

COMPLEMENT SIGN BIT OF EXPONENT. 

STORE IT. 

COUNT 24 MUL/23 DIV ITERATIONS. 

RETURN. 

IF POSITIVE EXP THEN NO OVFL. 



Wozpak ][, November 1979, pages 109-115. 



FLOATING POINT PACKAGE 



The mantissa-exponent, or 
widely used by computers t 
floating point representat 
memory to store than the n 
floating point arithmetic 
subroutine package in ROM, 
Maximum precision is retai 
such as 'divide by zero' a 
point number representatio 
as floating point BASIC. 



'floating point' numerical representation is 
o express values with a wide dynamic range. With 
ion, the number 7.5 x 10~22 requires no more 
umber 75 does. We have allowed for binary 
on the APPLE ] [ computer by providing a useful 
which performs the common arithmetic functions, 
ned by these routines and overflow conditions 
re trapped for the user. The 4-byte floating 
n is compatible with future APPLE products such 



A small amount of memory in Page Zero is dedicated to the floating point 
workspace, including the two floating-point accumulators, FP1 and FP2. 
After placing operands in these accumulators, the user calls subroutines in 
the ROM which perform the desired arithmetic operations, leaving results in 
FP1. Should an overflow condition occur, a jump to location $3F5 is 
executed, allowing a user routine to take appropriate action. 



FLOATING POINT REPRESENTATION 



HI 



LOW 



Exponent 
1. Mantissa 



Signed Mantissa 



The floating point mantissa is stored in two's complement representation 
with the sign at the most significant bit (MSB) position of the high-order 
mantissa byte. The mantissa provides 24 bits of precision, including sign, 
and can represent 24-bit integers precisely. Extending precision is simply 
a matter of adding bytes at the low order end of the mantissa. 

Except for magnitudes less than 2^-128 (which lose precision) mantissa are 
normalized by the floating point routines to retain maximum precision. 
That is, the numbers are adjusted so that the upper two high-order mantissa 
bits are unequal. 

HIGH-ORDER MANTISSA BYTE 
01.XXXXXX Positive mantissa. 
10.XXXXXX Negative mantissa. 
00.XXXXXX Unnormalized mantissa. 
ll.XXXXXX Exponent = -128. 

2. Exponent. 

The exponent is a binary scaling factor (power of two) which is applied to 
the mantissa. Ranging from -128 to +127, the exponent is stored in 
standard two's complement representation except for the sign bit which is 
complemented. This representation allows direct comparison of exponents, 
since they are stored in increasing numerical sequence. The most negative 
exponent, corresponding to the smallest magnitude, -128, is stored as $00 
($ means hexidecimal) and the most positive, +127, is stored as $FF (all 
ones) . 



EXPONENT 


STORED AS 


+127 


11111111 ($FF) 


+3 
+2 


10000011 ($83) 
10000010 ($82) 



+1 10000001 ($81) 

10000000 ($80) 

-1 01111111 ($7F) 

-2 01111110 ($7E) 

-3 01111101 ($7D) 

-128 00000000 ($00) 
The smallest magnitude which can be represented is 2^-150. 



HIGH LOW 

EXP MANTISSA 



The largest positive magnitude which can be represented is +2^128-1. 



$7F I | $7F I | $FF I | $FF 



EXP MANTISSA 



FLOATING POINT REPRESENTATION EXAMPLES 



DECIMAL 


HEX 


HEX 


NUMBER 


EXPONENT 


MANTISSA 


+ 3 


81 


60 


00 00 


+ 4 


82 


40 


00 00 


+ 5 


82 


50 


00 00 


+ 7 


82 


70 


00 00 


+12 


83 


60 


00 00 


+15 


83 


78 


00 00 


+17 


84 


44 


00 00 


+20 


84 


50 


00 00 


+60 


85 


78 


00 00 


- 3 


81 


A0 


00 00 


- 4 


81 


80 


00 00 


- 5 


82 


B0 


00 00 


- 7 


82 


90 


00 00 


-12 


83 


A0 


00 00 


-15 


83 


88 


00 00 


-17 


84 


BC 


00 00 


-20 


84 


B0 


00 00 


-60 


85 


88 


00 00 



FLOATING POINT SUBROUTINE DESCRIPTIONS 

FCOMPL subroutine (address $F4A4) 

Purpose: FCOMPL is used to negate floating point numbers. 

Entry: A normalized or unnormalized value is in FP1 (floating point 
accumulator 1) . 

Uses: NORM, RTLOG. 

Exit: The value in FP1 is negated and then normalized to retain precision. 
The 3-byte FP1 extension, E, may also be altered but FP2 and SIGN are not 
disturbed. The 6502 A-REG is altered and the X-REG is cleared. The Y-REG 
is not disturbed. 

Caution: Attempting to negate -2~128 will result in an overflow since 
+2~128 is not representable, and a jump to location $3F5 will be executed, 
with the following contents in FP1. 



FP1: I I I $80 



XI Ml 



Example: Prior to calling FCOMPL, FP1 contains +15. 



FP1: | $83 J | $78 J | | j | (+15) 

XI Ml 
After calling FCOMPL as a subroutine, FP1 contains -15. 



FP1: ] $83 I | $88 I | | | | (+15] 



XI Ml 

FADD subroutine (address $F46E) 

Purpose: To add two numbers in floating point form. 

Entry: The two addends are in FP1 and FP2 respectively. For maximum 



precision, both should be normalized. 

Uses: SWPALGN, ADD, NORM, RTLOG. 

Exit: The normalized sum is left in FP1. FP2 contains the addend of 
greatest magnitude. E is altered but sign is not. The A-REG is altered 
and the X-REG is cleared. The sum mantissa is truncated to 24 bits. 

Caution: Overflow may result if the sum is less that -2~128 or greater than 
+2^128-1. If so, a jump to location $3F5 is executed leaving in XI, and 
twice the proper sum in the mantissa Ml. The sign bit is left in the 
carry, for positive, 1 for negative. 



FP1: I I I X.YYY. 



XI Ml 
(For carry=0, true sum=+X.YYY x 2"128) 
Example: Prior to calling FADD, FP1 contains +12 and FP2 contains -5. 



FP1: | $83 | $60 I I I I (+121 



XI Ml 



FP2: | $82 I | $B0 I | | | | (-5) 



X2 M2 
After calling FADD, FP1 contains +7 (FP2 contains +12) 



FP1 | $82 | $70 | | I I (+7) 



XI Ml 

FSUB subroutine (address $F468) 

Purpose: To subtract two floating point numbers. 

Entry: The minuend is in FP1 and the subtrahend is in FP2. Both should be 
normalized to retain maximum precision prior to calling FSUB. 

Uses: FCOMPL, ALGNSWP, FADD, ADD, NORM, RTLOG. 



Exit: The normalized difference is in FP1 with the mantissa truncated to 24 
bits. FP2 holds either the minued or the negated subtrahend, whichever is 
of greater magnitude. E is altered but SIGN and SCR are not. the A-REG is 
altered and the X-REG is cleared. The Y-REG is not disturbed. 

Cautions: An exit to location S3F5 is taken if the result is less than 
-2^128 or greater than +2"128-1. or if the subtrahend is -2"128. 

Example: Prior to calling FSUB, FP1 contains +7 (minuend) and FP2 contains 
-5 (subtrahend) . 



FP1: I $82 I | $70 I 10 1 I I (+121 



XI Ml 



FP2: I $82 I | $B0 I j I I | (-5; 



X2 M2 
After calling FSUB, FP1 contains +12 and FP2 contains +7. 



FP1: I $83 I I $60 I 10 1 I I (+121 



XI Ml 

FMUL subroutine (address $F48C) 

Purpose: To multiply floating point numbers. 

Entry: The multiplicand and multiplier must reside in FP1 and FP2 
respectively. Both should be normalized prior to calling FMUL to retain 
maximum precision. 

Uses: MD1, MD2, RTL0G1, ADD, MDEND. 

Exit: The signed normalized floating point product is left in FP1. Ml is 
truncated to contain the 24 most significant mantissa bits (including 
sign). The absolute value of the multiplier mantissa (M2) is left in FP2. 
E, SIGN, and SCR are altered. The A- and X-REGs are altered and the Y-REG 
contains $FF upon exit. 

Cautions: An exit to location $3F5 is taken if the product is less than 
-2"128 or greater than +2"128-1. 

Notes: FMUL will run faster if the absolute value of the multiplier 



mantissa contains fewer 'l's than the absolute value of the multiplicand 
mantissa. 

Example: Prior to calling FMUL, FP1 contains +12 and FP2 contains -5. 



FP1: ] $83 | | $60 | | | | | (+12) 



XI Ml 



FP2: I $82 I | $B0 I I I I | (-5; 



X2 M2 
After calling FMUL, FP1 contains -60 and FP2 contains +5. 



FP1: | $85 | | $88 | | | | | (-60; 



XI Ml 



FP2: I $82 I I $50 i I I I I (+ 5] 



X2 M2 

FDIV subroutine (addr $F4B2) 

Purpose: To perform division of floating point numbers. 

Entry: The normalized dividend is in FP2 and the normalized divisor is in 
FP1. 

Exit: The signed normalized floating point quotient is left in FP1. The 
mantissa (Ml) is truncated to 24 bits. The 3-bit Ml extension (E) contains 
the absolute value of the divisor mantissa. MD2, SIGN, and SCR are 
altered. The A- and X-REGs are altered and the Y-REG is cleared. 

Uses: MD1, MD2, MDEND. 

Cautions: An exit to location $3F5 is taken if the quotient is less than 
-2"128 or greater than +2"128-1 

Notes: MD2 contains the remainder mantissa (equivalent to the MOD 
function). The remainder exponent is the same as the quotient exponent, or 
1 less if the dividend mantissa magnitude is less than the divisor mantissa 



magnitude. 

Example: Prior to calling FDIV, FP1 contains -60 (dividend), and FP2 
contains +12 (divisor) . 



FP1: ] $85 | | $80 | | | | | (-60) 



XI Ml 



FP2 I $83 I I $60 I 10 1 I | (+12! 



XI Ml 
After calling FMUL, FP1 contains -5 and M2 contains 0. 



FP1: | $82 | | $B0 | | | | | (-5) 



XI Ml 

FLOAT Subroutine (address $F451) 

Purpose: To convert integers to floating point representation. 

Entry: A signed (two's complement) 2-byte integer is stored in Ml 
(high-order byte) and Ml+1 (low-order byte) . Ml+2 must be cleared by user 
prior to entry. 

Uses: N0RM1. 

Exit: The normalized floating point equivalent is left in FP1. E, FP2, 
SIGN, and SCR are not disturbed. The A-REG contains a copy of the 
high-order mantissa byte upon exit but the X- and Y-REGs are not disturbed. 
The carry is cleared. 

Notes: To float a 1-byte integer, place it in Ml+1 and clear Ml as well as 
Ml+2 prior to calling FLOAT. 

FLOAT takes approximately 3 msec, lonqer to convert zero to floating point 
form than other arguments. The user may check for zero prior to calling 
FLOAT and increase throughput. 

* 

* LOW-ORDER INT. BYTE IN A-REG 

* HIGH-ORDER BYTE IN Y-REG 

* 

85 FA XFLOAT STA Ml+1 



84 F9 


STY 


Ml INIT MANT1 


A0 00 


LDY 


#$0 


84 FB 


STY 


Ml+2 


05 D9 


ORA 


Ml CHK BOTH 
BYTES FOR 


DO 03 


BNE 


TOFLOAT ZERO 


85 F8 


STA 


XI IF SO CLR XI 


60 


RTS 


AND RETURN 


4C 51 F4 TOFLOAT JMP 


FLOAT ELSE FLOAT 






INTEGER 



Example: Float +274 ($0112 hex) 
CALLING SEQUENCE 



A0 01 


LDY #$01 


HIGH-ORDER 
INTEGER BYTE 


A9 12 


LDA #$12 


LOW-ORDER 
INTEGER BYTE 


84 F9 


STY Ml 




85 FA 


STA Ml+1 




A9 00 


LDA #$00 




85 F8 


STA Ml+2 




20 51 F4 


JSR FLOAT 





Upon returning from FLOAT, FP1 contains the floating point representation 
of +274. 



FP1 I $88 I I $44 I $80 I I I (+274) 



XI Ml 



FIX subroutine (address $F640) 

Purpose: To extract the integer portion of a floating point number with 
truncation (ENTIER function). 

Entry: A floating point value is in FP1. It need not be normalized. 

Uses: RTAR. 

Exit: The two-byte signed two's complement representation of the integer 
portion is left in Ml (high-order byte) and Ml+1 (low-order byte) . The 
floating point values +24.63 and -61.2 are converted to the integers +24 
and -61 respectively. FP1 and E are altered but FP2, E, SIGN, and SCR are 
not. The A- and X-REGs are altered but the Y-REG is not. 

Example: The floating point value +274 is in FP1 prior to calling FIX. 



FP1: | $88 | $44 | $80 | | (+274) 



XI Ml 

After calling FIX, Ml (high-order byte) and Ml+1 (low-order byte) contain 
the integer representation of +274 ($0112) . 



FP1: I $8E | $01 I I $12 I 10 



XI Ml 
Note: FP1 contains an unnormalized representation of +274 upon exit. 

NORM Subroutine (address $F463) 

Purpose: To normalize the value in FP1, thus insuring maximum precision. 

Entry: A normalized or unnormalized value is in FP1. 

Exit: The value in FP1 is normalized. A zero mantissa will exit with X1=0 
(2 exponent). If the exponent on exit is -128 (X1=0) then the mantissa 
(Ml) is not necessarily normalized (with the two high-order mantissa bits 
unequal). E, FP2, SIGN, AND SCR are not distubed. The A-REG is disturbed 
but the X- and Y-REGs are not. The carry is set. 

Example: FP1 contains +12 in unnormalized form (as .0011 x 2 ). 



FP1: ] $86 | | $0C | ] | I I (+121 



xl Ml 
Upon exit from NORM, FP1 contains +12 in normalized form (as 1.1 x 2 ) 



FP1: I $83 I I $60 I 10 1 I I (+12) 



XI Ml 



N0RM1 subroutine (address $F455) 

Purpose: To normalize a floating point value in FP1 when it is known the 
exponent is not -128 (X1=0) upon entry. 

Entry: An unnormalized number is in FP1. The exponent byte should not be 
for normal use. 



Exit: The normalized value is in FP1. E, FP2, SIGN, and SCR are not not 
disturbed. The A-REG is altered but the X- and Y-REGs are not. 



ADD Subroutine (address $F425) 

Purpose: To add the two mantissas (Ml and M2) as 3-byte integers. 

Entry: Two mantissas are in Ml (through Ml+2) and M2 (through M2+2) . They 
should be aligned, that is with identical exponents, for use in the FADD 
and FSUB subroutines. 

Exit: the 24-bit integer sum is in Ml (high-order byte in Ml, low-order 
byte in Ml+2). FP2, XI, E, SIGN and SCR are not disturbed. The A-REG 
contains the high-order byte of the sum, the X-REG contains $FF and the 
Y-REG is not altered. The carry is the '25th' sum bit. 

Example: FP1 contains +5 and FP2 contains +7 prior to calling ADD. 



FP1: | $82 | $50 I I I I (+5) 



XI Ml 



FP2: | $82 | $70 | | | | (+7) 



Upon exit, Ml contains the overflow value for +12. Note that the sign bit 
is incorrect. This is taken care of with a call to the right shift 
routine. 



FP: | $82 | $C0 | | I I (+121 



ABSWAP Subroutine (address $F437) 

Purpose: To take the absolute value of FP1 and then swap FP1 with FP2. 
Note that two sequential calls to ABSWAP will take the absolute values of 
both FP1 and FP2 in preparation for a multiply or divide. 

Entry: FP1 and FP2 contain floating point values. 

Exit: The absolute value of the original FP1 contents are in FP2 and the 
original FP2 contents are in FP1. The least significant bit of SIGN is 
complemented if a negation takes place (if the original FP1 contents are 
negative) by means of an increment. SCR and E are used. The A-REG 



contains a copy of X2, the X-REG is cleared, and the Y-REG is not altered, 



RTAR Subroutine (address $F47D) 

Purpose: To shift Ml right one bit position while incrementing XI to 
compensate for scale. This is roughly the opposite of the NORM subroutine. 

Entry: A normalized or unnormalized floating point value is in FP1. 

Exit: The 6-byte field MANT1 and E is shifted right one bit arithmetically 
and XI is incremented by 1 to retain proper scale. The sign bit of MANT1 
(MSB of Ml) is unchanged. FP2, SIGN, and SCR are not disturbed. The A-REG 
contains the least significant byte of E (E+2), the X-REG is cleared, and 
the Y-REG is not disturbed. 

Caution: If XI increments of (overflow) then an exit to location $3F5 is 
taken, the A-REG contains the high-order MANT1 byte, Ml and XI is cleared. 
FP2, SIGN, SCR, and the X- and Y-REGs are not disturbed. 

Uses: RTLOG 

Example: Prior to calling RTAR, FP1 contains the normalized value -7. 



FP1 I $83 I I $A0 I 10 1 10 1 (-7) 



XI Ml 

After calling RTAR, FP1 contains the unnormalized value -7 (note that 
precision is lost off the low-order end of Ml) . 



FP1 $84 I $D0 10 I I (-7) 



XI Ml 
Note: Ml sign bit is unchanged. 

RTLOG subroutine (address $F480) 

Purpose: To shift the 6-byte field MANT1 and E one bit to the right (toward 
the least significant bit). The 6502 carry bit is shifted into the 
high-order Ml bit. This is useful in correcting binary sum overflows. 

Entry: A normalized or unnormalized floating point value is in FP1. The 
carry must be cleared or set by the user since it is shifted Into the sign 
bit of Ml. 

Exit: Same as RTAR except that the sign of Ml is not preserved (it is set 



to the value of the carry bit on entry) 

Caution: Same as RTAR. 

Example: Prior to calling RTLOG, FP1 contains the normalized value -12 and 
the carry is clear. 



FP1: I $83 I | $A0 I I I I | (-12; 



XI Ml 

After calling RTLOG, Ml is shifted one bit to the right and the sign bit is 
clear. XI is incremented by 1. 



FP1: | $84 | $50 | | | | (+20) 



XI Ml 

Note: The bit shifted off the end of MANT1 is rotated into the high-order 
bit of the 3-byte extension E. The 3-byte E field is also shifted one bit 
to the right. 



RTL0G1 subroutine (address $F484) 

Purpose: To shift MANT1 and E right one bit without adjusting XI. This is 
used by the multiply loop. The carry is shifted into the sign bit of 
MANT1. 

Entry: Ml and E contain a 6-byte unsigned field. E is the 3-byte low-order 
extension of MANT1. 

Exit: Same as RTLOG except that XI is not altered and an overflow exit 
cannot occur. 



MD2 subroutine (address $F4E2) 

Purpose: To clear the 3-byte MANT1 field for FMUL and FDIV, check for 
inital result exponent overflow (and underflow), and initialize the X-REG 
to $17 for loop counting. 

Entry: the X-REG is cleared by the user since it is placed in the 3 bytes 
of MANT1. The A-REG contains the result of an exponent addition (FMUL) or 
subtraction (FDIV). The carry and sign status bits should be set according 
to this addition or subtraction for overflow and underflow determination. 

Exit: The 3 bytes of Ml are cleared (or all set to the contents of the 
X-REG on Entry) and the Y-REG is loaded with $17. The sign bit of the 



A-REG is complemented and a copy of the A-REG is stored in XI. FP2, SIGN, 
SCR, and the X-REG are not disturbed. 

Uses: NORM. 

Caution: Exponent overflow results in an exit to location $3F5. Exponent 
underflow results in an early return from the calling subroutine (FDIV or 
FMUL) with a floating point zero in FP1. Because MD2 pops a return address 
off the stack, it may only be called by another subroutine. 

### 



